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Abstract 

A simple interlayer pair tunneling model is solved exactly. We find that in the 
normal state the spin-| particle and hole excitations are gapped. But the state is 
an unusual metal, characterized by novel fermionic, spin-zero and charge ib2e gapless 
excitations that exist about new type of Fermi surfaces. The model is consistent with a 
number of unusual properties of underdoped cuprates. Superconductivity is induced by 
an additional intra-layer pairing interaction which opens a gap in the charge spectrum. 
The symmetry of the order parameter is in general different from the symmetry of the 
single-electron (or hole) gap. The former can be a d-wave, while the the latter is more 
complex. 
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In underdoped bilayer cuprate superconductors there appears to be a gap, known as 
the spin gap, in the normal-state spectrum. For example, below about 150 K, magnetic 
susceptibility x ^'^^ NMR relaxation rate decrease rapidly Yet there is apparently no 
charge gap since the in-plane conductivity is metallic. Recent photoemission experiments 
further show that only a fraction of the the Fermi surface near (vr, 0) and symmetric points 
seems to be gapped 0. On continuity grounds, Fermi or Luttinger liquids are not expected 
to show such a behavior, suggesting that these systems may belong to a different universality 
class. Since the spin-gap has been observed only in bilayer systems, it has been suggested 
that interlayer interactions may be responsible 0. An interesting possibility is interlayer 
pair tunneling which was initially considered as a mechanism for enhancement 0. The 
underlying assumption is that single electrons can not hop coherently between copper-oxide 
layers so that second order hopping processes such as pair tunneling become important. 

Recently Chakravarty et al. have studied a simple model of pair tunneling between 
two layers and found considerable enhancement of T^. within a mean-field approximation 
0. Physics within each layer is described by a Bardeen, Cooper and Schrieffer (BCS) type 
reduced Hamiltonian, and tunneling is also confined to the reduced subspace of pairs with 
zero total momentum. The mean-field state does not have a spin gap. But, as shown by 
Anderson, when the tunneling interaction is treated exactly the normal-state x exhibits spin- 
gap behavior 0. In this paper we show that, although the model is too simple to provide 
a detailed description of the cuprates, it has a number interesting properties. We find that 
in the normal state spin-| electron and hole excitations acquire a gap. The Fermi surface 
appears partially gapped at finite T (experimental situation). An important property is 
that there are gapless spin-0, charge ±2e excitations which are fermionic and have their own 
Fermi surface, and which dominate the low-T physics. The system is thus metallic, but not a 
Fermi liquid in the usual sense. To induce superconductivity we add an intralayer attractive 
interaction which opens a gap in the charge spectrum. Therefore the symmetry of the charge 
gap (and hence the order parameter) is determined by the attractive interaction and is in 
general different from the symmetry of the one particle (i.e., spin-) gap. 

We consider the following Hamiltonian for two layers 

^ = IIe(k)4kaCika - 51 ^kk'clkt4-kiC*,-k'iCi,k'T - II^j(k)[4,kT4,-kiC2,-kiC2,kT + h.c.]. 

«k(T k,k' k 

(1) 

Here Q^ko- destroys an electron carrying a (two-dimensional) wavevector k, spin cr =|, | in 
layer i = 1,2. The first term describes an electron band of energy e(k) for each layer. The 
second term is an in-plane attractive interaction and the third term decribes tunneling. Like 
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the BCS Hamiltonian only pairs of zero total momentum are considered in the interaction 
terms so that tunneling matrix element TjCk) is diagonal in k. The assumption is that 
Tj(k) is 0(1). This is different from more traditional interactions such as Vkk' (or tunneling 
terms with nonzero momentum) which scale inversely as the volume of the system and hence 
does not contribute for k = k'. It is precisely the diagonal nature of Tj(k) that causes the 
enhancement of Tc and the spin gap. Although the problem can be solved for any Tj(k), for 
definiteness we will take 

Tj(k) = -^{cosk^ - cosfcy)^. (2) 

This is the same form used in the mean-field analysis [^]. Our focus will be on the universal 
features and as such we will take a simple cosine band: e(k) = —2t[cos k^a + cos kyo]. 

We first ignore V and solve the normal state problem. Then the Hamiltonian can be 
diagonalized exactly by diagonalizing each k subspace separately 0. For a given k, we 
need to consider four single-electron states: (l,k |),(1,— k |),(2,k t) and (2, — k J,). It 
is convenient to diagonalize = -f^k ~ /wA^k; where A^k is the number operator and /i is 
the chemical potential. Let ^k = e(k) — /U. There are 16 many-body states which can 
be grouped according to the total electron number A^k and z-component of total spin Sz 
which are both conserved. In the absence of tunneling, an eigenstate has energy A^k^k, with 
A'k = 0, 1,2,3,4. We represent such a state by \a b >, where a (6) stands for one of the 
four states in the first (second) layer: 0, |, J., ti- Tunneling only connects two of these: 
III > and |0 Ti>) with A'k = 2 and = 0. Diagonalizing the 2 by 2 matrix leads 
to the following two states: 

|2±>=^[| n 0> ± |0 Ti>], (3) 

with energies 2^k =p Tj(k), respectively. The remaining fourteen states are unaffected. Note 
that since A^k is conserved, there is no long-range order, i.e, there is no superconductivity in 
the absence of V. 

Ground State: The ground state of the full Hamiltonian is constructed by selecting the 
lowest energy state of Tik for each k. The relevant states are the (i) the four-electron state 
|4 >= I tl t|> of energy 4^k; (ii) the zero-electron state |0 >= |0 > of energy and 
(iii) the two-electron state |2+ > of energy 2(^k — Tj{k). For Tj = 0, states with different 
A^k cross at the Fermi level which is at ^k = 0. So one must pick the four-electron state for 
^k < and the zero-electron state for ^k > 0. 

For nonzero Tj{k), the four-electron state still has the lowest energy for < ~\Tj{'k) 
and the zero-electron state for ^k > \Tj(k). But in the region — |Tj(k) < ^k < \Tj(k), 
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the two-electron state \2+ > has the lowest energy. We will call this the two-electron 
region. Therefore the original Fermi surface has disappeared, except at points where the line 
Tj{k) — intersects the Fermi surface f^k = 0/ In its place, two new "Fermi" surfaces 
have appeared. Surface I is at = — |Tj(k) and separates the four-clcctron states from the 
two-electron states. Siu^facc II appears at ^k = |7j(k) and separates the two-electron states 
from the zero-electron states. As shown in Fig 1, the particle number A^k has discontinuities 
at the new Fermi surfaces whose topolgy depends on the form of Tj(k). In the quartic case, 
Tj(k) is largest near (0, ±7r) and (±7r, 0), but it vanishes along kx = ky, so that the original 
Fermi surface survives at four points. 

Excitations: In the absence of tunneling, the elementary excitations are the usual holes, 
created by picking a three-electron state in the four-electron region and particles, created 
by picking a one-electron state in the zero-electron region. These are gapless (energy |^k|), 
carry charge ±e and spin |. 

(i) For Tj(k) > 0, the hole and particle excitations described above exist for ^k < 

and ,^k > |^j(k)j respectively. Two additional spin-^ excitations exist in the two-electron 
region: a particle of energy ^k + Tj{k) and a hole of energy — ,^k + Tj{h). All four spin-^ 
excitations are gapped. Therefore, magnetic susceptibility x is zero at T = 0. The "spin gap" 
is thus a natural consequence of tunneling, will be observed in photoemission experiments. 
We define a gap function Aspj„(k) by the value of the excitation energy at the new Fermi 
surface. This equals ^Tj(k), exactly the distance between the old and new Fermi surfaces. 

(ii) Near the new Fermi surfaces there are additional gapless hole- and particle-type 
excitations. However, they have spin zero and charge ±2e. Near surface I, a particle is 
created by selecting the four-electron state in the two-electron region and a hole is created 
by selecting a two-electron state in the four-electron region. These have energy ±(2{k+7j(k)) 
which vanish at surface I (^k = ~\Tj{k)), and can not be created from two spin-| excitations. 
Similar spin-zero, charge-two particle and hole excitations also exist near Fermi surface II 
with energies ±(2^k — Tj{k)). 

In short, while there is a spin gap, there is no charge gap. The charge excitations can 
not be viewed as 'preformed' Cooper pairs since, as shown below, they behave like fermions, 
at least as far as thermodynamic properties are concerned. The partition function for the k 
subspace is given by 

Z(k) = Zo(k) + Zi(k) - (1 + e-^«'')^ + 2[cosh(/?rj(k)) - l]e-2/^«N (4) 

where Zq is the partition function for Tj — 0, and /? = 1/kT. Thermodynamic properties 
are largely determined by the behavior of the two factors Zq/Z and Zi/Z as a function of 
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T. At high temperatures, kT » Tj, Zi ^ and we recover the noninteracting result. But 
at low temperatures such that PTj » 1 and » 1, the two factors look like Fermi 
functions. Thus for < 

Zo/Z ~ /(2e + Tj), Z^/Z ~ (1 - /(2e + Tj)). (5) 

And for ^ > 

Zo/Z^{l - f{2C - Tj)), Zi/Z~/(2e - Tj). (6) 

Note that these functions have step discontinuities at the new Fermi surfaces, and the energy 
argument correspond to the charge excitations. Of course, far from the Fermi surfaces they 
lose their fermionic character. 

We will be interested in the case for which the noninteracting Fermi energy epo » tj, kT. 
As shown in Fig. 1, for kT << Tj, the new Fermi surfaces and the fermionic character of 
the charge excitations are evident in the particle number = which is given by 

Zo(k) Zi(k) 

With increasing T, the charge excitations break up and density of spin-| excitations increases, 
and at kT ~ Tj, the new Fermi surfaces have effectively disappeared so that nk looks like the 
Fermi distribution function for the noninteracting system. Therefore the gap Aspj„(k) will 
be seen only on those parts of the Fermi surface where Tj{k) is large compared to kT. For 
the quartic case, Tj(k) vanishes at four points along the (1,1) direction, and is quite small 
(because of the fourth power) except in the neighborhood of (ivr, 0) and (0, ±7r), where the 
gap is largest. Hence, for T > 0, parts of the Fermi surface along (1,1) and nearby directions 
will appear to be ungapped. For a fixed T, the gap decreases with doping, as the old Fermi 
surface moves away from the maximal points. For fixed density it effectively disappears for 
kT > tj. These results are in qualitative agreement with recent photoemission measurements 
i- 

The state is clearly not the usual Fermi liquid. This is true for arbitrarily small tj, and 
all densities and spatial dimensions (there is no quantum phase transition). Moreover, there 
is no broken symmetry. For kT << tj, charge fermions dominate and make the system a 
degenerate metal. With increasing T these give way to spin-| excitations, and eventually at 
kT ~ tj there is a crossover to an ordinary Fermi liquid. 

The susceptibility is determined by the spin-^ excitations, and is given by 



^ ^0 ^o(k) 
^ 2N^ Z{k) 



(e/35k + 1)2 



(8) 
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where /iq — eh/2mc. The quantity in the bracket is locahzed at the old Fermi surface = 0) 
and has width kT. For kT < Tj, its overlap with the factor Zq/Z is exponentially small. 
Thus if T'j(k) has no zeroes at the Fermi level (finite gap), then x vanishes exponentially 
with T. Suppose Tj(k) is a constant = 2Aspin- Then, for kT << Tj « epo, we find 
that x/Xo ~ {kT/Aspin)exp{—2Aspin/kT), where Xo = |poA*o the zero-temperature sus- 
ceptibility for the noninteracting system and po is the corresponding density of states at 
the Fermi level. For kT » tj, on the other hand, we recover the noninteracting result: 

X/XO ~ 1 - U^sp^n/kT)\ 

If the gap function has zeroes one expects a power law. The power depends on the 
symmetry of Aspj„(k). For the quartic case, Tj(k) (and hence the gap) vanishes at the 
Fermi points kg as Tj(k) oc {q^ — Qy)'^, where q = k — kg. Then we find that x vanishes 
as X ~ T^/^. More generally x ~ T^/p, if the spin gap vanishes with a power p. Fig 2. 
shows the temperature dependence of calculated numerically, for several values of the 
electron density n. We see that x deviates sharply from xo below a temperature Tsus and 
then decreases slowly as T^/^. The scale Tgus decreases with decreasing n. 

For T — > charge-fermions dominate. For example, specific heat per site C{T) = T^{T) 
vanishes linearly with T. Fig. 3 shows the temperature dependence of 7(T)/7o, where 70 is 
the coefficient for the noninteracting system. We find that at low T, 7(T)/7o ^ 1/4 + 07"^/^. 
The constant term 1/4 arises because, for charge fermions, the spin-degeneracy is reduced 
by a factor of 2 and their Fermi energy is incresed by a factor of 2, to leading order mtj/eFo- 
The T^/^ is again the contribution from the spin-| excitations. As shown in Fig. 3, ^{T) 
goes through broad maximum at T = Tsputi and at large T approaches the nonointeracting 
value 7o. Thus the Wilson ratio (oc CjTx) is temperature dependent and becomes infinite 
at T = 0. The scale Tg^^i is in general much smaller than Tgus- 

Superconductivity : To induce superconductivity we now add the in-plane interaction 
Vk,k'- It is sufficient to treat this term by mean- field approximation since it is restricted 
to the reduced subspace of zero-momentum pairs. But Tj is again treated exactly. Let 
fok =< c^kT^.-kj. > be the order-parameter which we take to be real and independent of the 
layer index. Then the mean-field pairing Hamiltonian is 

Hsuper,MF = - ^1 4^^-^ + Q -kiQ^kT , (9) 

i,k 

where = ^kk'^k'- The full Hamiltonian can again be diagonahzed for each k sepa- 
rately. The total is still a good quantum number for each k, but A^k is no longer conserved. 
States with Sz are not affected by tunnehng, and therefore are the same as in the BCS 
theory. In particular, the Sz — ±1/2 subspaces spht into two manifolds of energy 2^k ± -Rk, 
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where it^k = (Ck + ^D' ■ 

The lowest-energy eigenstate belongs to the same — subspace as before, the subspace 

consisting of the three states |0 >, |2+ > and |4 >. These are now coupled by the pairing 
Hamiltonian leading to the following 3 by 3 matrix: 

-V2Ak 

-V2Ak 2a-Tj(k) -^Ak 

-^Ak 4Ck 



The eigenvalues are given by 



Am = S^k 



2 1 

cos[-(27rm + 0)], 



(10) 



where m = 1,2,3, and tan(/) = - [^-c'^Y /c, d ^ A^l + AAl + |rj(k) and c = 



Tj{k) 2AI + |Tj(k) - Ail 



a. Charge Gap: For Tj = we have the usual BCS problem. The spin-| particles and holes 
acquire a gap Ak at the original Fermi surface. For Tj ^ 0, these are already gapped in the 
normal state. The main source of the condensation energy is a gap that opens up in the 
charge fermion spectrum. For example, |0 > and |2+ > are degenerate with energy 0. at 
the (new) right Fermi surface (^ = Tj/2). This is split by pairing, giving a charge gap which 
has the form AchargeO^) = Ak/i(2^^^), where fi{x) is even in x. For Ak small compared to 
2Tj(k), we find that AchargeO^) ~ 2-\/2Ak. The superconducting Tc is enhanced relative to 
the BCS value in part because AchargeO^) is generally larger than Ak. 

b. Svmmetrv of the Order Parameter: The order parameter 6k =< cik-fCi-kj. > is an odd 
function of Ak since a change of sign Ak — > — Ak is undone by a gauge transformation 
Ckf — > — Ckt which changes the sign of 6k- At T = it has the form 



6k = Al 



(4e-A(k))-^-A-^(k) 
l + 2A2(A(k)-2 + (4e-A(k))- 



(11) 



where A(k) is the lowest eigenvalue of the 3 by 3 matrix and is even in Ak. The symmetry of 
the order parameter (and of the charge gap) is thus determined by the symmetry of the in- 
plane gap function Ak, which is determined by our choice of Vk,k'- In particular it could have 
d^2_j^2 symmetry. On the other hand, the gap in the spin-| electron and hole excitations — 
the analog of the BCS gap — does not have the same symmetry since it now depends both on 
Tj(k) and Ak. In other words, experiments that probe the order parameter (e.g., Josephson 
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tunneling) could find a d-wave symmetry, but experiments in which single electrons are 
added or removed (e.g., photoemission) would see a more complex gap structure. 

A number of other results can be obtained independently of the origin or details of Vk,k'- 
As in the BCS model, Vk,k' is presumed to be appreciable near the original Fermi surface (at 

= 0), i.e., for ^k,^k' < i^o, where uq is the analog of the Debye energy. Consequently Ak 
decreases with distance from ^k = 0. But the charge gap scales with the value of Ak at the 
new Fermi surface which is at a distance |Tj(k) = Aspj„(k). It follows that AchargeiM) is small 
where Agpinik) is large (e.g., near ivr, 0). This implies that close to half filling (underdoped 
regime) charge gaps are small, and increases with doping (decreasing n). Therefore, in this 
region, (1) Tc also increases with doping, as observed. (2) For fixed T << Tc, density of 
charge excitations decreases with doping. Consquently, electronic contribution to specific 
heat, thermal conductivity etc would decrease with increasing doping. (3) Since Tc increases 
and Tsus decreases with doping, the temperature range over which the spin gap is seen 
decreases with increasing doping and eventually disappears. 

In conclusion, the interlayer pair tunneling model provides a simple example of the break 
down of Landau theory and can account for a number of unusual properties of the underdoped 
cuprates. Since the effect is caused by a k-diagonal tunneling interaction one may ask 
whether such a restriction is physically justified. For example, a more general tunneling 
Hamiltonian which results from second order hopping processes would have the form 

= IZ^^(k,q)[cI,kT^l,-k+qiC2,-k+qiC2,kT + h.c], (12) 

k,q 

where q is the total momentum of the pair. Now, extensitivity of free energy requires that 
Sk,q^j(k, q) scales with the volume V, so that T7(kq) oc V~^. But Hamiltonian (1) is 
obitained if only q = terms are kept. This is not thermodynamically significant unless one 
postulates that Tj(k, q = 0) is 0(1), whereas terms with q 7^ scale as and hence are 
similar to other interaction processes. Note that this implies the existence of something akin 
to a condensate which favors zero-momentum pair tunneling. Phenomenologically, if the di- 
agonal terms are responsible for the spin-gap state, then thermodynamic considerations alone 
guarantee that they are fundamentally different from other interaction processes. It then 
makes sense to start with diagonal terms only, and treat the remainder as fluctuatuations. 
Generalization to many layers is, in principle, straightforward. Indeed Baskaran has shown 
that, if V is ignored, the model can be mapped onto a solvable quantum spin-chain model 
along the c-axis for each k, provided the Hamiltonian is restricted to the A^k = even (i.e., 
charge) subspace only 0. Then it is not clear what happens to the spin gap. Unfortunately, 
when V is included the spin-chain model is no longer solvable. 
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Figure Captions 

Fig. 1. Particle density N]^ along the (tt, 0) direction. The intermediate step with N^k = 2 
is the two-electron region. At low T = O.ltj {ks = I) there arc two Fermi functions repre- 
senting charge ±2e excitations. But at T = 0.5tj, there is only one Function corresponding 
to the nontinter acting FS (which is where the three lines intersect). Bandwidth — 40ij. 

Fig 2. (a) Magnetic susceptibility x/Xo as a function of scaled temperature t/tj, for four 
values of electron density n. Note that characteristic temperature decreases with decreaing 
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n. (b) At low T, X vanishes as T^/^, as shown in the log-log plot. The slope of the hne is 
1/4. 

Fig 3. (a) Scaled specific heat ccfiicient ^(T)/^o — C{T)/Co{T) as a function of of 
scaled temperature T/tj. The broad maximum appears at a temperature which is much 
lower than temperature at which susceptibility drops (see Fig 2a) and may not be seen 
if superconductivity intervenes, (b) The same quantity on a log-log plate showing that it 
approaches a constant (=1/4 -|- corrections) as T ^ 0. 
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